Intermittent dynamics and 1//^ noise in single cardiac muscle cells 
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Fluctuations in the spontaneous beating activity of isolated cardiac cells were studied over a 
timescale of six decades. The beat dynamics of single cardiac cells were heterogeneous and intermit- 
tent. The interbeat intervals (IBIs) were power- law distributed in a long-time regime. Furthermore, 
for long timescales up to the experimental window, the autocorrelation of IBIs exhibits a scaling 
behavior of l//^-noise type. These observations suggest that 1//^ noise is an intrinsic characteristic 
of spontaneous activity of single cardiac cells. 
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Recently, it is rapidly recognized that noise is ubiq- 
uitous in cellular biology. Owing to the development of 
single- molecule experiments on various kinds of proteins, 
their stochastic motion has been extensively studied [1]. 
It has also been observed that several cellular behaviors, 
including gene expression [2] and motile activity [s*], ex- 
hibit large noise. The properties and the significance of 
such cellular noise are currently of great interest in both 
physics and biology. 

It is also well known that a variety of signals from mul- 
ticellular organisms exhibit characteristic fluctuations [1] . 
A typical example of such fluctuations is that observed in 
human heart-beat rate [5,l6l]. It is now widely accepted 
that the power spectrum of the human heart beat rate 
exhibits 1//^ dependence {jS 1) in the low- frequency 
regime. However, in spite of intensive studies in the last 
two decades [7] , the mechanism oil/ noise generation 
still remains unclear. 

In the preset Letter, in order to reveal the properties of 
the cellular noise in cardiac cells and to understand their 
relevance to the fluctuations at the multicellular level, we 
study the fluctuations in the activity of cardiac muscle 
cells at the single-cell level. Although there have been 
several studies reporting that the tissue culture of the 
cardiac muscle cells also exhibit scale-invariant fluctua- 
tions [1, @, d^, dll , it is unknown whether 1//^ noise is 
an intrinsic property of a single isolated cardiac cell. 

In order to clarify this point, we studied the sponta- 
neous activity of single isolated cardiac muscle cells in 
a wide range of time scales, ranging from subseconds to 
several days. It is then found that the spontaneous beat- 
ing activity of single isolated cells exhibits strong hetero- 
geneity and intermittency. The probability distribution 
of the interbeat interval (IBI) exhibits a power-law de- 
pendence in the long-time regime. Furthermore, a long- 
range correlation analysis, termed detrended fluctuation 
analysis (DFA) [lH, indicates that the time series of IBI 
fluctuations exhibit scale invariance over a timescales of 
three decades and are characterized as 1//^ noise. This 
experimental finding suggests that the 1//^ noise is char- 
acteristic of an isolated single cardiac cell. 

Spontaneous beating activity was examined for single 



isolated ventricular cells derived from neonatal rats. The 
cells were plated on a petri dish, which is coated with col- 
lagen, at a low density (30 cells/mm^) [13]. At this den- 
sity, most of the cells have no contact with other cells (see 
Fig. [1]). The cardiac muscle cells, which were identified 
by their morphology and ability to display spontaneous 
contraction, were isolated from each other during the ex- 
periment. The culture dish was set in a customized in- 
cubation chamber placed on an inverted microscope, and 
the cells were cultivated at a constant condition (37° C, 
5% CO2, saturated humidity). Phase-contrast images of 
cells were obtained via a charge-coupled device camera 
at 30 frame/s and were stored on a hard disk for the 
following analysis. 

Contractions of a cardiac muscle cell cause a temporal 
change in the brightness of the phase-contrast images. 
Thus, the timings of the spontaneous contractions of the 
cell were determined from the time series of the bright- 
ness of a pixel in the cell The inset in Fig. [T] displays 
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FIG. 1: Phase-contrast image of isolated cells in culture. The 
cells whose spontaneous activity was studied are marked by 
solid squares. The scale bar represents 100 /xm. Inset: A part 
of the time series of 8-bit brightness for the cell indicated by 
the white arrow (top). Raster plot constructed from the time 
series of brightness (bottom) . The timings of the spontaneous 
contractions are marked by vertical lines. 
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FIG. 2: Typical raster plot of a single cell. Each vertical line represents the timing of a single contraction. The numbers 
associated with this plot indicate the time elapsed since the onset of measurement. 



a typical example of a time series of brightness of a single 
cell, presented with a raster plot (timings of contractions) 
determined from the original time series. By employing 
the above mentioned experimental setup, the profiles of 
spontaneous activities of 11 cells from three different ex- 
periments have been analyzed. For each cell, data were 
obtained for several hours (mean 32.5 h, ranging from 6 
to 55 h). 

Figure [2] displays a typical example of the raster plot 
of a single cell. As observed in this plot, the timings of 
contractions are heterogeneous. We found that there are 
several typical patterns in this trace. 

First, steady contractions with some fluctuations were 
frequently observed, as observed in the portion around 
110000 s in Fig. [21 This pattern is represented in Fig. [3^ 
and is termed pattern A. It is also found that contrac- 
tions often become intermittent. Bursts of contractions 
are separated by quiescent periods. This pattern is ob- 
served around 40000 s in Fig. [21 for instance, and is shown 
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FIG. 3: Magnified raster plots of several typical temporal pat- 
terns; (a) Pattern A shows steady contractions with rather 
small fluctuations and (b) pattern B depicts an intermittent 
bursting pattern. The durations of the active and quiescent 
phases vary stochastically, (c) Pattern C shows a regular 
bursting pattern. Active and quiescent phases occurs alter- 
nately with a regular periodicity. 



at a higher resolution in Fig. [3)3. This pattern is termed 
B. In this pattern, the duration of quiescent periods oc- 
casionally reaches several hours. We also notice that in 
some cases, the bursting activity becomes quite regular. 
This pattern is represented in Fig. [3t. In this case, qui- 
escent periods occur with high regularity, although the 
periods of alternation vary slowly. This pattern is termed 
C. 

Although we can categorize the patterns of the contrac- 
tion as mentioned above, this categorization seems to be 
rather subjective. In the entire time series for the con- 
tractions, a pattern is gradually replaced with another 
pattern, and subsequent patterns do not have distinct 
boundaries. It might be more appropriate to define the 
patterns in terms of the duration of the quiescent peri- 
ods: small or no quiescent periods (pattern A), quiescent 
periods with durations ranging from several seconds to 
several minutes (pattern B), regular quiescent periods 
(pattern C). 

In order to characterize the statistics of the IBIs, we 
calculated their distribution. Figure [4] displays the dis- 
tribution of IBIs for the same cell as that in Fig. [2] on the 
logarithmic scale. It is found that a peak appears around 
0.5 s, which indicates the presence of a refractory period 
shorter than 0.5 s. Furthermore, for larger timescales, 
a scaling behavior is observed over a timescale of two 
decades. 

The exponent of scaling seems to depend on the pat- 
terns of activity of each cell. For cells that exhibit a large 
fraction of pattern B (intermittent), the exponent is close 
to 2. For cells exhibiting a large fraction of pattern A 
(regular), the scaling behavior becomes less distinct and 
the exponent increases. Moreover, for cells that exhibit 
a large fraction of pattern C (regular burst), the second 
peak appears around several tens of seconds, which is in- 
terpreted as the modal length of the quiescent periods. 
The mean value of the exponent averaged over the all 
cells was —3.1 ± 1.2. 

Next, in order to study the two-point statistics of the 
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FIG. 4: Frequency of IBIs in a logarithmic plot. A fit in the 
region [10"°•^ 10^-^] yields the exponent -2.23 =b 0.02. 



IBI time series, we performed a long-term memory analy- 
sis (the DFA) [12']. This method enables us to investigate 
the autocorrelation of a time series by removing the ef- 
fect of nonstationarity. First, letting represent the i-th 
IBI, the IBI time series are integrated as Xi = Xl}=i '^j- 
Subsequently, the entire time series is divided into m sub- 
sections of size n; in each of these subsections, the linear 
local trend is calculated by means of the least-square fit. 
After the local trend is subtracted in each subsection, the 
variance of the residual time series is calculated. There- 
fore, the quantity to be analyzed here is represented as 

m (j+l)n 2 

where represents the local trend in the j-th subsection. 
If the time series has scale invar iance, this quantity scales 
with n as 

V{n) (X n^. (2) 

The exponent a obtained in this manner is related to the 
scaling exponent /3 of the power spectrum of the time 
series Xi as a = /3 + 1. 

Figure [5] displays a typical example of V{n)^ which is 
calculated for the same cell as that chosen to plot Fig. [4j 
plotted against n on the logarithmic scale. Two scaling 
regions are prominent in this plot. One is the short- 
time region (n < 10^), where the exponent is a :^ 1. 
The other is the long-time region (n > 10^), where the 
exponent is a 2. Therefore, the fluctuations in the IBI 
time series are uncorrelated for short time scales (shorter 
than 10^ beats). For time scales larger than 10^ beats, 
it has been found that the IBI time series has a power- 
law correlation, and this correlation can be characterized 
as 1//^ noise. This power-law correlation seems to be 



FIG. 5: (Color online) Typical result of the DFA for a single 
cell. Solid circles represent V{n) calculated for the original 
data obtained in the experiment. The solid line is a fit in the 
region [10^,10^"^], and its slope yields the scaling exponent 
a = 1.79 ± 0.03. Solid squares represent V{n) calculated for 
the surrogate (randomly shuffled) data, and its slope yields 
the exponent close to unity (a = 0.96 ±0.02). Inset represents 
the frequency of the exponents in the long-term scaling region. 



sustained until at least 10^ beats, although this limit is 
set by the period of the experiment (55 h). 

For comparison, we prepared a randomly shuffled IBI 
time series as surrogate data. As expected, the surrogate 
data yields V{n) with a single scaling exponent, close to 
unity for all the time scales. From Fig. O it is apparent 
that the slope of the original data is distinct from that 
of the surrogate data. 

In order to examine the cell-to-cell variance of the 
abovementioned tendency, the frequency of a is presented 
in the inset of Fig. [5l The value of the exponent is deter- 
mined by fitting log V{n) to a log n + 6 in the scaling re- 
gion n > 10^. From this figure, it is found that the value 
of a is around 2. The mean value of a was 1.97 ± 0.23. 
This suggests that the 1 //^-noise-type long-term corre- 
lation is a generic property observed in the activity of 
single cardiac muscle cells. 

The abovementioned experimental findings suggest 
that intermittent, power-law distributed, and power-law 
correlated IBIs are the inherent properties of a single car- 
diac muscle cell. Here, let us summarize earlier studies 
concerning the fluctuations of IBIs in a cardiac tissue. 
Kucera et al. reported that the beat rates exhibit a 1//^ 
spectrum in a dense network of cardiac muscle cells, al- 
though it was not clear in their report whether it is an 
intrinsic property of single cardiac muscle cells or a phe- 
nomenon that arises through interactions among cells Q . 
Yoneyama et al. performed the DFA for isolated cardiac 
cells and reported that the exponent a/2 is close to 0.5 
in the short-time regime (n < 10^), which is consistent 
with our result [10]. 



4 



Soen et al. studied the beatrate fluctuations in single 
isolated cardiac cells in terms of the Allan Factor (AF) 
measure [9]. In their report, they reported that the IBIs 
are uncorrelated in the time scale t < 10^ x (mean IBI), 
which is consistent with our result; however, IBIs exhibit 
strong anti-correlation for larger time scales. This ob- 
servation corresponds to a 3 for the DFA for a long 
time scale, which differs from our result. One possible 
reason for this discrepancy might be the difference in the 
quantity analyzed. The essential difference between the 
DFA and the AF measure is that the former deals with 
the beat number as a measure of time, while the latter 
deals directly with the physical time. When applied to 
a highly heterogeneous time series, the DFA is preferred 
to AF measure. The difference between both measures is 
evident when applied to the surrogate data. As observed 
in Fig. [5l the DFA yields nice scaling with a :^ 1 for the 
all time scales, while the AF provides signiflcantly less 
clear scaling for the surrogate data (see Ref. [9|). 

In conclusion, it has been found that a single cardiac 
muscle cell exhibits heterogeneous and intermittent dy- 
namics. The IBIs obey a power-law distribution, and 
the time series of the IBIs are power-law correlated and 
exhibit 1//^ noise in the long-time regime. This correla- 
tion is sustained for at least several tens of hours. These 
properties of single cardiac cells might affect the heart- 
beat dynamics in vivo. 

The mechanism behind these singular properties at the 
single-cell level is an open problem. The fact that 1//^ 
noise is stably observed implies the presence of a robust 
mechanism for these intermittent fluctuations. With re- 
gard to this problem, it has been reported that membrane 
currents in several types of cells exhibit a power-law cor- 
relation [15| . In the investigation of the mechanism of the 
power-law correlation in single-cell fluctuations, cardiac 
muscle cells have a great advantage because one can in- 
vestigate the IBI time series from their contractile motion 
in a noninvasive manner, and this fact facilitates parallel 
and extremely long-term measurements. This advantage 
will enable us to clarify the fundamental mechanism be- 
hind power-law correlations in the cellular activity at the 
single-cell level. 
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